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Abstract. This paper suggests a new way to compute the path integral for simple 
quantum mechanical systems. The new algorithm originated from previous research in 
►^ ' string theory. However, its essential simplicity is best illustrated in the case of a free 

^^ , non relativistic particle, discussed here, and can be appreciated by most students taking 

l^ ' an introductory course in Quantum Mechanics. Indeed, the emphasis is on the role 

played by the entire family of classical trajectories in terms of which the path integral 
is computed exactly using a functional representation of the Dirac delta-distribution. 
On ' We argue that the new algorithm leads to a deeper insight into the connection between 

>^^ ■ classical and quantum systems, especially those encountered in high energy physics. 
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1. Introduction 

It is often argued tliat a formulation of quantum mechanics in terms of path integrals is 
too advanced to he within the scope of most undergraduate courses. On the other hand, 
a great deal of physics is done, nowadays, using Feynman's path integral method, with 
applications ranging from gauge theories in high energy physics to solid state physics, 
and statistical mechanics. As a result, the amount of pedagogical literature dedicated to 
a discussion of the path integral approach to quantum mechanics is increasing steadily 

iSiii- 

Quantum mechanics is probably the most challenging paradigm of physics that a 

student will ever come across, and the pedagogical effectiveness of Feynman's method, 

beautifully expounded by the author himself in his original work 0, lies in its appeal 

to intuition while addressing the most fundamental principles of the theory. Indeed, it 

seems to us that one of the major achievements of Feynman's formulation of quantum 

mechanics was to restore the particle trajectory concept at the quantum level. Our 

purpose, then, is to suggests a specific algorithm that emphasizes the privileged role 

that classical trajectories play in constructing the "sum over histories" envisaged by 

Feynman in order to describe the quantum evolution of a particle. 

This paper is an outgrowth of an earlier investigation in string theory where the 

computational method described here was first applied J^, |, ||. However, in order 

to illustrate the method in the simplest context, we apply it to the familiar case of a 

free, non relativistic particle. We show that the path integral can be computed exactly 

by summing over the whole family of trajectories which are solutions of the classical 

equations of motion. This is in contrast to the conventional approach to the path 

integral in which Newtonian dynamics singles out a unique classical trajectory which is 

then perturbed by quantum fluctuations. The main ideas and techniques are discussed 

in the case of a free particle with the bare minimum of formal apparatus which should 

be within grasp of most undergraduate physics majors taking a first course in quantum 

mechanics. 

Hopefully, the applicability of our approach to higher dimensional objects will justify the 

claim, made in the abstract, that our method applies to "simple quantum mechanical 

systems" besides the point particle case considered in this paper. 

In order to place our work in the right perspective, we begin Section ^ with a brief 

overview of the discretization procedure that forms the basis of the "sum over histories" 

approach to Quantum Mechanics. This will give us the opportunity to single out the 

precise point of departure of our own approach from the conventional one based on the 

"Lagrangian path integral" . 

In Section ^ we construct the propagation kernel and make contact with the diffusion 

equation of a free, non relativistic particle. 

In [Appendix A| we outline a possible extension of our definition of path integral to the 

case of a particle in an external potential. 

Finally, [Appendix B| is a compendium of three basic formulae which are especially 
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relevant to the discussion in the text. 

2. Some Remarks about Path Integration 

In this section we review briefly the deflnition of phase space path integral. The main 
purpose here is twofold: flrst, we wish to present a self-contained exposition for those 
readers who may not be acquainted with the idea of "sum over histories" proposed by 
Feynman|^, |T^; second, we wish to isolate those fundamental properties of Quantum 



Mechanics, encoded in the path integral approach, which are relevant to the subsequent 
discussion of our own computational method . 

One such fundamental property is the distinction between the square of the absolute 
value of the wave function, or probability density, and the wave function itself, or 
probability amplitude. With that distinction in mind, the peculiarity of the quantum 
mechanical world is that, in contrast to the classical rules of conditional probability, 
one computes probability amplitudes for the paths, and then sum the amplitudes; when 
amplitudes are superimposed, and then squared, an interference pattern is induced in 
the probability density. The famed "particle-wave duality" of the quantum world stems 
essentially from that new computational rule. According to Dirac and Feynman |[11 



the basic difference between classical and quantum mechanics is that the former selects, 
through the stationary action principle, a single trajectory connecting the initial and 
flnal particle position, while the latter assumes that a particle "moves simultaneously 
along all possible trajectories" connecting the initial and flnal end points. Of course, 
not all trajectories are equal, even if they correspond to wave functions having the 
same absolute value: the action S[y(t)] of a particular path connecting the two 
fixed end points determines the phase of the propagation amplitude associated with 
that path. The phase factor corresponding to any individual path can be written in 
the form exp {i S [y(t) , p(t)] /h}, and a very effective way to describe the quantum 
mechanical behavior of a particle is through the sum over histories proposed by 
Feynman. Symbolically, the amplitude to evolve from Xq to x is 

Aixo,x)= Yl ^[yit).pit)] , (1) 

paths from xq to x 

where 

<|.[j;(t),p(i)] = /C-exp|^5[j/(t),p(t)]| . (2) 

A sum of this type, deflned over a space of functions, i.e., y(t), and p(t), is a 
functional integral, and the basic problem of the Feynman formulation of quantum 
mechanics is to determine what that formal sum means. 

In order to construct the sum over all paths in phase space, one may follow the evolution 
of the system in a total time interval [t ,t'] by dividing such interval in N+1 subintervals, 
with end points labeled by ti, ^2, ■ ■ ■ , ti , . . . , t^+i so that 

t =to 

t = tN+l 
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For convenience, one may also choose each of the above intervals to be of equal length 
e, 



(3) 
(4) 



ti+i — ti — e , 
in such a way that 

t^ = t + ke , k = 0,...,N + 1 . 

Next, one has to determine an equation for the partial amplitude 

^ ( Qi+i ! ti+i ',qi,ti) = { Qi+i , ti+i I 5'j , tj ) (5) 

that the system evolves from a position with coordinate qi at time tj to a position with 
coordinate g^+i at time tj+i. Using the bracket notation, the amplitude for the first 
subinterval takes the form 

(g,e|g,0) = (g|exp(-i/Je) \q) 

= 6{q-q)-ie{q\H\q) + 0{t^) 

dp ir>(a-a) f ^P i„(a-a) : . rr ( .. ^ + ^ 



27r . 2tt '^ 



dPM-^- 

2n 


-Q) 


dp 
2vr^"P 


{• 



1 — itH p 



q + q 



2 



0{e^) 



p\ 



q)-eH [ p, 



q + q 



Generalizing the above expression to a generic subinterval, we find 



U+i,ti+i\qi,ti) = I — exp<j? 



dpi 



Pi 



i+1 - qi) - H [ Pi, [ti+i - ti 



(6) 



(7) 



Approximating the path from q at time t to q' at time t' with a finer subdivision of the 
total time interval, and in view of the fact that there are no predetermined conditions 
on the coordinate variable at intermediate times, the total amplitude can be written 
as the product of the amphtudes associated with each subinterval, integrated over all 
possible intermediate positions: 

^ ^ ^ l,Af+l 

A{q' ,t' ■,q,t) = ^lip^ dqi dq2... dqN Y\ ("?»' ^»5 g^-i, ti-i) • (8) 

Combining the two equations above, we are led to the final step of the discretization 
procedure 



A(q't':q,t)= lim 



dqi / dq2 
dpi f dp2 



dqN- 
dpN+i 



271 J 2n J J 2tx 

.g4p2(<?2-gi)-(t2-tl)J/(p2,^4^)] .^__. 
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'1,N „ \ /1,N+1 
N-t+oo 



/1,N „ \ /1,N+1 



/hN \ /1,N+1 



' ' 3 



^^i'"^^ 



„ 1k-1k-l rr/'„ Qk+lk-l ' 
Pk f, _f, — --n-\pk, o 



(tk—tk-l) 



The expression in the last equahty represents the discrete version of the path integral 
in phase space 

J [Vq] j [Vp] exp {i[pq-H {p, q)] dt} . (9) 

which is the starting point for computing the sum over histories for the majority of the 
integrable systems encountered in the literature. Equation (||) is also the starting point 
of our own approach. However, in terms of the discretization procedure outlined above, 
a further simple elaboration of the phase space path integral leads to a computational 
method which seems mathematically more efficient and physically more enlightening, 
at least for some complex systems encountered in high energy physics. 
Consider the first term in the expression (^, namely 

pqdt = pdq . (10) 



This term corresponds to the discrete sum 

1,JV+1 

J^ Pk {qk - Qk-i) (11) 

k 

which we can rewrite as follows: 

1,N+1 

Y^ Pkiqk- Qk-i ) = 

k 

= X^ PkQk - 5Z Pk^i'-'^ 

k k 

1,N+1 1,N+1 1,N+1 1,N+1 

= Y PkQk - Y PkQk^i - 5Z Pk-i9k-i + "Y Pk-iQk-i 

K K K K 

1,N+1 1,N+1 

= ^ {pkQk - Pk-iQk-i) - Y (Pk- Pk-i) Qk-i 
k k 

1,N 

= PN+iqN+i-poqo-Y{pk-Pk-i)qk-i ■ (12) 
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The first two terms in the above expression represent a boundary contribution to the 
path integral, whereas the remaining sum corresponds to the following integral 

^{pk-Pk-i)qk-i^ qdp= qpdt . (13) 

Next, in what follows we make use of the following equality 

dqk-i exp {i {pk - Pk^i) qk-i} = S {pk - Pk-i) (14) 



which gives the well known representation of the Dirac delta function as Fourier 
transform of the imaginary exponential. Thus, the net result of the above rearrangement 
of terms is encoded in the following correspondence 

y"[Pg]exp|^y"rftgp|=5[p] =^ jim^AT (t,-) H^ (f^^^) • (15) 

which we take as a definition of the "functional Dirac delta distribution" . 
In the next section, we apply this mathematical rearrangement of terms in the phase 
space path integral to compute the propagation kernel of a non relativistic particle. As 
we shall see, the physical payoff is a novel interpretation of the sum over histories in 
the sense that it underscores the special role played by the entire family of trajectories 
which are solutions of the classical equations of motion. 

3. Non-Relativistic Particle 

3.1. Path Integral and Propagation Kernel 

The classical dynamics of a non-relativistic point-particle with mass m is encoded into 
the Lagrange function, or in its corresponding Hamiltonian 

L(y(t),^(t);t)=im^2 (^g) 

H{m,p{t);t) = ^p' . (17) 

/^From here follow the classical equations of motion: 

dp ^ ^ , ^ 

-j- = => p = cost. = q (18) 

^ = ^ H = cost. = E (19) 

dt 

^^ = Q ^ yit) = -Qt + xo , (20) 

dt m 

where we have taken into account the boundary conditions 

j/(0) = 5o , y{T) = ^ => q = m^^^ (21) 



so that 



y{t) = ^^j^t + Xo . (22) 
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Furthermore, for later reference, we recall that the classical action 

Sc\X ^ , Xq ; T ) = -m / dt 



2 io \ T 

m , ^ ^ i9 

is a solution of the Jacobi equation 

AC 1 

^ + ;^V^ci. ■ V^ci. = . (23) 

Equivalently, one can follow the particle evolution by determining the propagation 
amplitude K{x — Xq; T) that a particle propagates from an initial position Xq to a 
final position x. As we have shown in the previous section, the amplitude is given by 
a formal sum over all phase space paths connecting Xq to a; in a total time T, each 
path carrying a weight given by the phase factor exp {i S[path.]/h). We have also 
reviewed the standard method that gives meaning to the sum over histories: it goes 
through a discretization procedure which consists in subdividing the total time lapse 
T into a number of infinitesimal intervals, thereby approximating the smooth phase 
space trajectory followed by a classical particle with a succession of "jagged" paths. It 
seems worth emphasizing at this point, that in the above implementation of the path 
integral, the classical trajectory followed by a particle is uniquely specified by Newtonian 
mechanics. Furthermore, we hasten to say that the sum over histories, so defined, is 
a sum over trajectories which are nowhere differentiable. As a matter of fact, a post 
modern interpretation of Feynman's discussion is that the quantum mechanical path 
of a particle is inherently fracta^ The gist of the argument is that, when a particle 
is more and more precisely located in space, its trajectory becomes more and more 
erratic as a consequence of Heisenberg's principle. In other words, it is the addition of 
quantum fluctuations around the classical trajectory that gives meaning to the idea that 
"a particle moves along all possible paths" connecting the initial and final configurations. 
Be that as it may, the discretization procedure turns the functional integral into an 
infinite product of ordinary integrals [^ . For some elementary systems, the integration 
over the momenta can be carried out, albeit with some efforts, and the final result for 
the amplitude is the "Lagrangian path integral" : 

/•* i r^ 

Kix,xo;T)= I [Vyit)]exp- L[y{t),m;t]dt . (24) 

For more complex systems, such as the relativistic extended systems encountered in 
contemporary high energy physics, the expression (P^, may not tell the full story, and 
we have found it advantageous to start with a sum over histories in phase space. Thus, 
in order to illustrate our computational method, we start directly from the canonical 

"*" This interpretation, clearly implied in the book by Feynman and HibbspO|, was revisited many years 
ago by Abbott and WiseQ, and reviewed by the authors in the larger, but mostly unexplored context 
of the microscopic structure of quantum spacetimett p[. 
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phase space path integral for a free, non-relativistic particle 



K{x,Xo;T' 
where 



A^ 



Hoip) 



X 



X, 



[Vy{t)][Vp{t)]exp- 



2m 



dt 



pit)-m-Hoip) 



(25) 



(26) 



and A^ is a normalization constant to be fixed later on. 

The focal point of our approach is the recognition that the path integral automatically 
assigns a special role to the whole family of trajectories which are solutions of the 
classical equation of motion. The precise meaning of the above statement is illustrated 



by implementing the computational steps following Eq. 
term in the action as 



Thus, we write the first 



dtp-y 



dt 



d_ 
Jt 



[py\ 



dty-p 



'0 ^0 "'' ^0 

Accordingly, the path integral reads 
K{x, Xo;T) = 



d{p-y) - dtyp 



Xo 



(27) 



N 



X 



Xo 



[Vy{t)][Vp{t)]exp{- 



X 



d{p-y) 



Xn 



dty{t)-p{t) 



dtHoip) 



and we note that, while the first term in the above expression is the integral of a total 
differential, and therefore independent of the path connecting the two points {xq,x), 
the second term depends linearly on the spatial trajectory y{t). A closer look at this 
term shows that it is a (functional) Dirac-delta distribution represented as a "Fourier 
integral" over the functions y{t) 



X 



[Vy{t)]exp{-- I dty{t)-p{t) 



6 



dp 
lit 



{2t 



'Xo k "-JO 

The Dirac-delta in equation (^) is non- vanishing only when its argument is zero, i.e., 
when the momentum p solves the classical equation of motion 

-^ = => p{t) = const. = q . (29) 

This is our first, and central result which applies equally well, "mutatis mutandis", to 
relativistic point-particles and relativistic extended objects: once the spatial trajectories 
y{t) are integrated out, the resulting path integral is non-vanishing only when the 
surviving integration variable, namely the three-momentum vector along the trajectory, 
is constrained to satisfy equation (^). This information is encoded in the following 
expression for the propagation amplitude 

K{x, Xo;T) = 



N [Dp{t)]S 



dp 



exp 



X 



nk^^-^ 



exp <^ -- / dtHo{p) 



(30) 
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where only the restricted family of constant (i.e., time independent) p-trajectories 
contributes to the path integral. Note that the value of the three-momentum along a 
classical trajectory is a fixed, even though arbitrary, number. Thus, on the mathematical 
side, the pay-off of our procedure is that we have traded the original path integral, i.e., 
a functional integral, with a single ordinary integral over the constant, but numerically 
arbitrary, components of the three-momentum. 

On the physical side, our method underscores a conceptual difference between the 
mechanism of propagation envisaged here and the conventional one. By this we mean 
that in the conventional interpretation of the path integral, the sum over histories is 
obtained by summing over all possible quantum fluctuations around a single classical 
trajectory which is uniquely defined because both extremal position and momentum are 
preassigned. In contrast, in our interpretation of sum over histories, only the extremal 
coordinates of the particle are precisely specified by the boundary conditions, whereas 
the corresponding value of the momentum at the end points is constant but arbitrary, in 
consistency with the uncertainty principle. Since the momentum is a vector, it follows 
that all paths in configuration space contribute to the evolution of the wave function. 
In the above sense, it seems to us that the complementarity between particle and wave 
behavior and the role of the uncertainty principle are especially evident in our approach. 
To summarize our discussion so far, "the sum over histories" of a free particle between 
two fixed end points, can be translated into the integration over all possible values of 
the linear momentum. Put briefiy, 

[Dp{t)] 6 [dp/dt] (...) = sum over the classical momenta 

= sum over constant momentum trajectories 

d^q (...) = ordinary momentum integral . 



Next, if the momentum is constant, then the Hamiltonian is constant as well, and we 
can write 



d'^q exp 









^1/2 



exp 



2Kr 



X — Xq 



where we have used the identity 

rX 

dy = X- xo 



,(31) 



(32) 



Xn 



and the formula for the multidimensional Gaussian integral listed in Appendix [Appendix 

O. As a result, we obtain 



K{x,Xo;T] 






im 
2Kf 



X — Xq 



(33) 



Finally, the normalization constant N can be fixed by the condition that, in the limit 
of a vanishing lapse of time, the particle is bound to be in its initial position. In other 
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words, 



lim K( X — Xq] T 



N 



2TT!f 



3/2 



10 



(34) 



Hence, 



K{x-Xo;T)= (^- 



m \3/2 
;) exp 



2nT' 



X 



Xo 



.2i7ihT. 
m \3/2 

^2iTrhTy 
where the phase factor is just the classical action measured in h units 



^— — j exp[iSd.{x-xo]T)/h] 



Sci.ix- Xo;T] 



ra 
~2 



dt 



X — Xq 

T 



m 

2r' 



X 



Xo 



(35) 



(36) 



The above expression for the propagation kernel is a well known result which tests the 
consistency of our approach: as Feynman and Hibbs demonstrated, for any lagrangian 
quadratic in the position and velocity variables, the corresponding propagation amplitude 
is given by a pre-f actor, which is a function of the evolution parameter lapse, multiplied 
by the exponential of "i" times the classical action in h units. 

3.2. Diffusion Equation and Propagation Kernel 

In the previous subsection, we have deduced the form of the propagation kernel (l35|), by 
evaluating the sum over histories as a Gaussian integral. With this result in hands, one 
can make contact with the more familiar formulation of quantum mechanics by showing 
that the propagation kernel K satisfies a diffusion equation of the Schrodinger type. To 
see this, first we calculate 



^A-(a--«„;T)=A'/<f',(-.«j)exp 


p- 1 dy 

" JXo 


exp j 


' iT 
2hm 


and 






/\^K{ x-Xo;T) = N d^q ( -^— ^ ) exp 

J \ h^ J 


^q- dy 


exp< 


iT 
2hm* 


Comparing the two equations above, we conclude that 






h"^ d 

— A,ir( x-Xo;T)= ih—K 


'{x-Xo;T) 







q q 



qq 



(37) 



(3^ 



(39) 
Finally, recalling the relation between the amplitude K and the wave function, namely 

'ip{x, t) = d^XoK{x - Xq; t)'ilj{xo, 0) , (40) 

we also conclude that ip{x ,t) must satisfy the time-dependent Schrodinger equation. 
At this point, it seems pedagogically instructive to reverse the procedure, and show that 
the propagation kernel can be determined by solving the diffusion equation ( p9|) which 
we now assume as given. To this end, we make the following ansatz, 

K{x-Xo;T) = F{T) exp{iM{x-Xo;T)/h) , (41) 
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in terms of two trial functions F{T) and M{x — Xq;T). The overall normalization 
constant may be determined by the same boundary condition (|34D. 
In order to determine the form of the trial functions, we demand that the tentative 
expression ( ^T]) satisfies equation (^). Thus, one finds 

^A.M(.-..;T)^-^ (42, 

and 

_L(VM).(VM)^-^ . ,«, 



Note that equation (|43|) is just the classical Jacobi equation p3|). Thus, without further 
calculations we can identify the phase of the kernel with the classical action: 

M{x-Xo;T) = S,i,{x-Xo;T) = —\x-Xo\^ . (44) 

Next, in order to determine the dependence of the kernel pre-f actor on T, we make use 
of equation (|^). Therefore, we first apply the Laplacian operator to SdX x — Xq; T ): 

—SdXx-Xo;T) = -ix~xo), (45) 

^^5ei.(5-5o;T) = -5., (46) 

S^i-^Sdix-xo;T) = '^6^ = ^ . (47) 

Then, we substitute the result into equation (0), which now takes the form 

dF 3 ^ dF 3dT 

— = F => — = . (48) 

dT 2T F 2 T ^ ^ 

Finally, integrating the last equation, we obtain 

F(r) = const. X r-^/2 . (49) 



Substituting all of the above in the original ansatz (0), leads to the following expression 
for the propagation kernel, 

K{x-Xo;T) = j;^exp{iS,,Xx-Xo;T)/h) . (50) 

Except for the normalization constant which is fixed by the initial condition, this is the 
same expression obtained from the path integral. 

4. Concluding Remarks 

As argued in Ref. [Q, the conceptual foundations of quantum mechanics are still the 
subject of debate, and thus "it is important to introduce undergraduate physics majors 
to non-standard developments of quantum mechanics" in the hope that exposure to 
such developments, early in the curriculum, might motivate some students to focus on 
the physical principles underlying the theory. In this paper we have taken up that 
challenging task by discussing a novel approach to the computation of the sum over 
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histories in the path integral. 

On the mathematical side, we have suggested a phase space formulation of the Feynman 
sum over paths in which the spatial trajectory y{t) plays the role of a Lagrange multiplier 
enforcing the classical equation of motion ( p!8[ ) at the quantum level. Imposing such a 



constraint, with the aid of a Dirac delta-distribution, reduces the functional integral to 
an ordinary integral over the arbitrary values of the initial momentum. In the simple 
case of a free, non relativistic particle, the use of the phase space path integral, or 
the Lagrangian formulation of it, is largely a matter of choice, in the sense that one 
arrives at the same expression of the propagation kernel. In the case of relativistic 
extended systems, the use of the phase space path integral along the same lines discussed 
here, seems definitely advantageous, if not mandatory. On the physical side, however, 
our formulation of the path integral reflects a new mechanism of propagation which, 
in our view, sheds some new light on the connection between classical and quantum 
systems. The element of novelty can be summarized thus: it is well known that classical 
mechanics, which is a deterministic theory, demands that one definite value of the 
initial position and momentum be assigned a priori, and thereafter selects a single 
classical trajectory that satisfies the equation of motion. On the other hand, in quantum 
mechanics one cannot fully specify the value of the momentum of a precisely localized 
particle. The price to pay for specifying the particle position at Xq when t = 0, is that p 
is completely arbitrary, and therefore one has to sum over all possible values of the three- 
momentum. This brings into the game the whole family of trajectories which satisfy the 
classical equations of motion. With all such channels of propagation open, a quantum 
particle turns to a "wave-like" propagation mode, exploring all of them by "evolving 
simultaneously" along all classical trajectories corresponding to all possible values and 
orientations of the momentum, up to the final point x where the particle is detected. 
Thus, as stated in most textbooks of Quantum Mechanics, but never fully clarified, a 
quantum particle starts as a pointlike object precisely located at a point, propagates as a 
wave, only to reappear as a dot on a screen, at the point of detection. We have exploited 
this dual, wavelike, description of the family of classical trajectories by connecting the 
quantum evolution of a particle with the Jacobi formulation of Classical Mechanics, i.e., 
by deriving the diffusion equation from the propagation kernel, and vice versa. The 
main purpose of the whole exercise was to test the consistency of our approach, and 
its ability to reproduce some standard results obtained by different methods. Equally 
important, however, the objective was to clarify, with the familiar example of a free 
particle, the connection between the classical theory and the central idea upon which 
the path integral formulation is based, namely, that "trajectories" still play a role in 
quantum mechanics and that all paths, with their intrinsic randomness, contribute to the 
evolution of the wave function which is governed by the propagation kernel. Having said 
that, we may add that our mathematical analysis and its physical interpretation extend 
well beyond the case of a free, non relativistic particle. Indeed, the very possibility of 
formulating a quantum mechanics of closed strings and other extended objects^, points 
to the effectiveness of our approach to the path integral for studying new fundamental 
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physics. By an interesting feedback process, we have found that many new formal and 
physical properties of our approach to strings can be applied equally well to the case 
of a point particle which is now regarded as the object of lowest dimensionality in a 
hierarchy of geometric objects to which the same dynamical principles seem to apply. 
iFrom the present introduction, the interested reader may eventually progress to a more 
advanced discussion involving relativistic extended systems, such as strings, membranes 
and p-branes, with a change in notation rather than substance. 

Appendix A. Particle in an External Potential: Short Distance Limit 



The results obtained in the previous sections for a free particle are, of course, exact. 
However, the same approach can be used to derive the short distance approximation for 
the particle's kernel in an arbitrary potential V{x). 

If |a; — :ro| is much smaller than the range over which V{x) varies significantly, then we 
can Taylor expand V{x) in the neighborhood of Xq'- 



Vix) ^Vo-F-ix-xo)~ 
VV{x)] ^ ^ . 

In this case, the path integral reads 
K{x-Xo;T) = N [ [Vy{t)]pPit)] x 



(A.l) 
(A.2) 



Xo 



X exp < — 
n 



X 



d{p-y\ 



T 



dty{t) 



dp{t) 
dt 



dtiH-F-{y-y,) 



..3) 



'Xo Jo "'' ^0 

where H = Hq + Vq. Integration over the particle trajectory gives a "Dirac-delta" whose 



argument yields the classical equation of motion 

dt 
which leads, in turn, to the expression for the classical momentum 

p{t) = Ft + q . 

Thus, the integration over the momentum trajectory can be defined as follows 
dp 



[Vpit)] 6 



dt 



-F 



...)^ d'q [Vpit)] 5 p{t)-Ft-q 



(A.4) 
(A.5) 

(A.6) 



and the resulting path integral takes the form 



K{x-x^-T\ 



( 



\27rff 



/ m 
V 2iTThT 




X 



h \ 2hm 



3/2 iX 

exp — I Vo 



2rp2 



F X 



F'T 
6m 



X 
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X exp 



irriT ( ^ ^ T^ -l 

X — xqA — - — F 



2hT \ 2hm 



(A.7) 



where the norinahzation constant has been chosen in such a way as to reproduce equation 
(^) once the potential is switched off. 

Appendix B. Three Basic Formulae 

One needs only three basic formulae in order to derive the main results reported in this 
paper: one is the extension of the Gaussian integral to a Z)-component vector variable, 
and the other two are representations for the Dirac delta-function. 





Generalized Gaussian Integral 








f ( \ ■ ■ ■ \ 

/ rfX exp -- X' Aij X^ +B'Xi = 
J \ 2 J 








= const. X ( detA )^/^ exp (b ■ A"^ 


■B/2) 


(B.l) 


where 










A = ( Aij ) and B = {Bi) 




(B.2) 



Dirac-delta representations 



a) Functional Fourier Transform: 












5[X(A)] = (const.) X f[VP{X) 


exp - 
V 


-,/< 


iAP(A) 


X(A)) 


(B.3) 


b) Gaussian Representation: 












S (x) = lim — exp (— x' 


7e) 








(B.4) 



[1] Some useful annotated references can be found in the Resource Letter by Gutzwiller M C 1998 

Am. J. Phys. 66 304 
[2] Holstein B R 1998 Am. J. Phys 66 583: this article discusses a number of techniques employed 

in the path integral method, and contains some useful references to the particle propagator in a 

linear potential 
[3] English L Q and Winters R R 1997 Am. J. Phys. 65 390 
[4] Cohen S M 1997 Am. J. Phys. 66 537 
[5] Dowrick N J 1997 Eur. J. Phys. 18 75 
[6] Feynman R P 1948 Rev. Mod. Phys. 20 367 
[7] Ansoldi S, AuriUa A and Spallucci E 1997 Phys. Rev. D56 2352 
[8] Ansoldi S, Aurilia A and Spallucci E 1999 Chaos, Solitons and Fractals "Superstrings, M, F, S 

. . . Theory" , 10 197 
[9] Ansoldi S, AuriUa A and Spallucci E 1996 Phys. Rev. D53 870 

Auriha A, Spallucci E and Vanzetta I 1994 Phys. Rev. D50 6490 

AuriUa A, Smailagic A and Spallucci E 1993 Phys. Rev. D47 2536 



Particle Propagator and Path Integral Derivation 15 

[10] Feynman R P and Hibbs A R 1965 Quantum Mechanics and Path Integrals (New York: Mc 
Graw-HiU) 

[11] Dirac P A M 1958 The Principles of Quantum Mechanics, (London: Oxford University Press) 

Feynman R P, reference M. This paper, as well as Dirac's original work can be found reprinted 
in : Selected Papers on Quantum Electrdynamics (New York: Schwinger J ed., Dover, 1958) 

[12] Abbott L F and Wise M 1981 Am. J. Phys. 49 37; see also, 

Nottale L 1992 Fractal Spacetime and Microphysics (World Scientific) 

[13] The path integral approach is discussed in many textbooks besides that by Feynman and Hibbs, 
reference [M. See, for instance 

Felsager B 1981 Geometry, Particles and Fields (Odense: Odense University Press) 
Ramond P 1981 Field Theory, A Modern Primer (Reading, Massachusetts: Benjamin/Cum- 
mings Publishers) 

Ryder L H 1985 Quantum Field Theory (Cambridge: Cambridge University Press) 
Schulman L S 1981 Techniques and Applications of Path Integration (New York: Wiley) 
Swanson M 1992 Path Integrals and Quantum Processes (New York: Academic Press) 



